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Abstract 



We analyse an extension of the minimal supersymmetric standard model including the dominant trilin- 
ear and bilinear R-parity violating contributions. We take the trilinear terms from the superpotential and 
the bilinear terms from the superpotential as well as the scalar potential. We compute the neutrino masses 
induced by those couplings and determine the allowed ranges of the R-parity violating parameters that are 
consistent with the latest SNO results, atmospheric data and the Chooz constraint. We also estimate the 
effective mass for neutrinoless double beta decay in such scenarios. 
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I. INTRODUCTION 



It is firmly established that non-zero neutrino masses can provide a solution to the observed solar 
[1,3 and atmospheric neutrino deficits ||. This requires two large mixing angles (very recent SNO data 
[1| favour the large mixing angle solution, LMA, for the solar anomaly) and two hierarchical neutrino 



mass squared differences (Am s 2 ohr <C Ar 



We assume that the three mass eigenstates are given by 



to, (i = 1,2,3) and we parametrise the rotation matrix from neutrino flavour (/) to mass (i) eigenstates as: 
Vfi — -K23(023)-Ri3(#i3)-Ri2(0i2)- The MNS matrix (analogous to the CKM mixing matrix for the quark sector) 
neglecting CP phases is given by [[|, 



Vfi = 



C12C13 
-C23S12 - C12S13S23 
S23S12 - C12C23S13 



C13S12 Sl3 
C12C23 - S12S13S23 C13S23 
-C12S23 - C23S12S13 C13C23 



(1) 



where & 



cos 6ij and s; 



sinOij. For the solar anomaly, which we take to be a consequence of f e -f p 



oscillation, the relevant mass squared difference is Amf 2 



solar' 



while for the atmospheric case, the 



oscillation being between and u T , the relevant mass squared difference is Amf 3 « A 



'23 



= Ar 



the inclusion of the SNO data, the MSW-LMA oscillation is the most favoured solution with Ar 
(2.5 - 19.0) x 1CT 5 eV 2 and sin 2 29 12 = 0.61 - 0.95 



After 



olar 



= (2 - 5) x 10" 3 eV 2 with sin 2 26> 23 = 0.88 - 1.0 



The SuperK atmospheric neutrino data suggest 
]. The Chooz || and Palo Verde long baseline 



reactor experiments bound sin 2 #13 0.04. Concerning the absolute masses, the recent claim for the evidence 
of neutrinoless double beta decay {Qv(3(3) by the Moscow-Heidelberg Collaboration || , although not yet firmly 
established, constrains the (ee)-element of the neutrino mass matrix in the flavour basis to lie between 0.1 and 
0.5 eV. Tritium /3-decay requires m v & 2.2 eV ||. Cosmological analysis from the recent 2dF Galaxy Redshift 
Survey constrains X)i TO i ^ 2.2 eV Q. 

Neutrino masses can be generated in the R-parity violating (Ifc p ) Supersymmetric Standard Model, where 
R p is defined as (_i) 3B + i + 2S 11 Here B,L and S are the baryon number, lepton number and spin of 
a particle, respectively. Strict phenomenological bounds on B and/or L violation exist in the literature |L3| ]. 
In this paper we assume that B violating couplings are absent and generate neutrino Majorana masses via 
two units of L violation. For this purpose, we allow both bilinear (/jj) and trilinear (A, A') interactions in the 
superpotential as well as the bilinear soft terms (£?,•) 0, 



W = [i J H u Lj + X -\ JKl LjL K El + \' J ^LjQ p D c q + h™H u Q p U c q , 
where the vector Lj = {Hd,L{) with J : 4..1, and the soft supersymmetry-breaking potential is 



(2) 



+A u p s H u Q p U^ + A j p s LjQ p D c s + -A JKl LjL K Ef + h.c. 



(3) 



As has been extensively discussed in refs |14| - |17| (and references therein), field redefinitions of the Hd,Li 
fields correspond to basis changes in Lj space and consequently the Lagrangian parameters will be altered. 
We use the basis-independent parameters constructed in Jl5|,[l6| and write the neutrino mass matrix in terms 
of 6^, 5 Z B , 5^ k , 5y q , which in the basis in which the sneutrino vevs are zero correspond to mVImIi ^VI-^Ij \jk, 
X' ipq , respectively. 

Calculations of neutrino masses in the context of lfc p theories have focused on tree-level contributions from 
the bilinear /ii-parameter or loop contributions from the trilinear A and A' couplings |18|. Recently a detailed 
phenomenological analysis has been done for a model including both A, A' and the superpotential bilinear 
parameter /ij (l^j2^] . In the basis- independent approach, a phenomenological analysis has been done including 
only the purely bilinear contributions, 8^ from the superpotential and 8 % fr om the soft Lagrangian [0. The 
possibility that loops involving 8^ were important was also discussed in ]17| . 



1 We follow closely the conventions given in Pq-[f7[ and we refer the reader to these papers for more details, 
symbols used in Eqs. (Q) and (^) have their usual significance. 
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Our purpose in this paper is to perform for the first time an analysis in a model which includes both 
the superpotential and soft bilinear parameters 6* and 5 B along with the superpotential trilinear couplings 

6^ , in the mass insertion approximation. A similar analysis can be performed including the trilinears S x p , q . 
In our analysis we take all fy p parameters to be real. Thus we update the analysis in jl7j by including 
the trilinear loop contributions alongside the bilinear tree and loop terms. We interface the neutrino mass 
matrix constructed out of the Ifip parameters with the constraints on the solar and atmospheric mass squared 
splittings, the corresponding mixing angles, and the Chooz constraint on #13. We then examine the nature of 
the mass spectrum and we check the consistency with cosmological data. We also obtain the effective mass 
for neutrinoless double beta decay in our scenario. Furthermore, we observe the impact of adding the trilinear 
contributions to the numerical results obtained in [[17f. 



II. PARAMETRIZATION OF THE MASS MATRIX 



The Ifip model we are considering generates a single neutrino mass at tree level proportional to 8^6 j^. 

There are additional loop corrections to the mass matrix when non- vanishing S B , S¥ (and 5 l f, q ) are included, 
leading to more than one non-zero mass eigenvalue. This enables us to fit the data on mass squared splittings 
and mixing angles. The relevant types of loops we consider in the mass insertion approximation are: 

• the well-known loops involving the trilinear I/L p couplings A or A' at the neutrino vertices I and II in 
Fig. 1 (with lepton/slepton or quark/squark as propagators), which give contributions proportional to 
S x Sx (or SySy); 

• the Grossman-Haber diagrams [0|, in which there are gauge couplings at the neutrino vertices while 
there are two types of $ p interactions contributing to the AL = 2 Majorana mass in the diagram of 
Fig. 0. The first kind have $ p couplings located at positions III + IV (slepton-Higgs mixing) with 
contributions proportional to SbSb- In the second type, the J]L P interactions are located at positions V 
+ IV (neutrino-neutralino and slepton-Higgs mixing) with contributions proportional to S^Sb', 

• the diagram of Fig. ^|, where two units of L violation come from positions V (neutrino-neutralino mixing) 
and II (A or A' vertex). The contribution to the neutrino mass is proportional to S p Sx (or S^Sx 1 )- 

The last two types of loops have frequently been overlooked in analytic estimates of neutrino masses. 
Exact formulae for these diagrams can be found in |l6| ] . Our formulae are robust (see later) and good enough 
for an order-of-magnitude estimate. 



Ill IV 



V 



I II 

FIG. 1. The conventional loops (Ijtp at I + II) and the Grossman-Haber loops ($ p at III + IV or V + IV) 
contributing to the neutrino mass. 



Setting all unknown sparticle masses equal to Ms — 100 GeV, and neglecting the mixing angles among 
the sparticles, we have a neutrino mass matrix of the form 



2 



v : 



ii 



FIG. 2. Charged loops with one gauge and a trilinear Yukawa coupling. The Ijtp interactions are at V and II. 



where 



lu]ij = M s 



» ^ cos/3 



cosV B B 



K2 



K3 



m e ,m ek 5T k ${ kn + 3 £ m d „m dfc ^ fc ^" 



L k 



(4) 



64tt 2 ' 



K2 



8tt 2 M s 



16tt 2 \/2' 



(5) 



We have included for completeness the contributions arising from the Sy terms which we put to zero in our 



numerical analysis. The simplest case to consider with a common 5 l u 



8b and 5 



ink 



5\ does 



not work as it gives only two non-zero masses and an eigenvector of the form (l/y/2, — 1/V2, 0) which cannot 
accomodate two large mixing angles for 9x2 and $23- 

In our numerical analysis we take : 5^, 5 B , 8™ k = 5\, (5™ fc = 0, for i = 1, 2, 3, i.e., seven independent 
parameters. Thus, the neutrino mass matrix elements will be given by, 



mn = M s 
mi2 = M s 
m 2 2 = M s 
m 13 = M s 
m 2 3 = M s 
m 3 3 = M s 



U2 



Kl 



;(^) 2 



cos 2 (3 

Kl Xl X2 _,_ K l 

cos j3 



Kl (S 2 B ) 2 + 2^S 2 ,5 2 



cos 2 /3 



cos/3 



4- 2K 3 m T Sj l S\ + K2fn T S x , 

+ K 3 rn T (5^5\ + S 2 8 x ) + K 2 mlSl, 
+ 2n 3 m T 6 2 5x + K 2 m 2 .Sl, (6) 



where we have employed the hierarchy of the charged fermion masses to keep only the dominant terms. 
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III. RESULTS AND CONCLUSIONS 



We have performed a general scan of parameter space made up by the seven parameters that appear in 
the mass matrix allowing tree- level contributions to either dominate over the loop corrections, to be on the 
same order as these, or to be much smaller than the loop terms. The fitted values of the couplings, obtained by 
using the atmospheric and Chooz data together with the preferred solar MSW-LMA solution from SNO, are 
presented in table I. Here we have taken cos (3—1 primarily to ensure a comparison with results in jljj ], which 
too employed the same value of cos f3, on the same footing. We stress that although this choice of cos (3 leads 
to an unacceptably low tan /3, we still use this value for illustration and effective comparison with previous 
results, noting at the same time that an order-of-magnitude estimate of the couplings is not sensitive to this 
choice. We also show in table I the allowed ranges of J2i m «i an d m cS = J2i Ki m » = J2i \ Vei\ 2 mi (since we 
have assumed the V matrix to be real) . This last quantity is the effective mass relevant for neutrinoless double 
beta decay (the neutrino masses induced by Tfl p interactions are Majorana type). 



Couplings 


Min 


Max 


Sx 


-2.0 x 10 -4 


2.0 x 10~ 4 


*i 


-6.8 x 1GT V 


6.8 x 10"' 




-8.4 x 10"' 


8.4 x 10 _Y 




-8.4 x 1Q-' 


8.4 x 10"'' 


si 


-2.7 x 10~ b 


2.7 x 10~ b 


5% 


-3.0 x 10~ b 


3.0 x 10~ 5 


5% 


-3.0 x 10" b 


3.0 x 10~ b 


m cfi (eV) 


1.9 x 10 -3 


6.7 x 10~ 2 


^rm(eV) 


4.9 x 


0.2 



TABLE I. Fitted range of the couplings satisfying MSW-LMA, SuperK and Chooz simultaneously. The fit is per- 
formed taking cos (3 = 1. 



In Fig. H we present the allowed region in the \Sb\ — (^s) 2 versus 1^1 = ^/X^(^i) 2 pl ane f° r the 

combined fit taking cos/3 = 1. We present our results for both S\ ^ and 5\ = 0. It can be clearly seen that 
the allowed region increases when we admit non-zero values of 5\. This is mainly due to the presence of the 
5fj,6\ terms in the mass matrix (originating from Fig. |^) which can take either sign and thus can accomodate 
a larger region of parameter space. 

The resulting fit strongly prefers a hierarchical mass pattern in our scenario, although our analysis cannot 
make a distinction between the inverted and the normal hierarchy. The inclusion of a non-zero S\ to the set of 
non-vanishing and 5 l B allows us to have all three non-vanishing neutrino masses while if only the latter two 
vectors are non-zero there are only two non-zero eigenvalues [jlTj . The maximum value of m e ff we predict can 
be tested in the next generation of neutrinoless double beta decay experiments. We have also checked that the 
cosmological bound restricting the sum of eigenvalues to be at most ~ 2 eV is fulfilled. 

To conclude, we outline the novel features of our analysis: (i) the trilinear ]ft, p interactions along with both 
types (originating from supcrpotcntial and soft terms) of bilinear Iji v terms have been employed for the first 
time to analyse the neutrino mass matrix; (ii) the seven parameter neutrino mass matrix has been interfaced 
with the latest SNO results favouring MSW-LMA solution along with the SuperK atmospheric and Chooz data; 
(iii) by including a non-vanishing 5\ we have updated the analysis in |l7| ] and the numerical impact of adding 
this term has been highlighted in Fig. ^; (iv) the mass spectrum in this kind of scenario implies a hierarchical 
pattern, although a distinction between the normal and inverted pictures cannot be made; (v) our prediction of 
m e ff can be tested in the next generation of neutrinoless double beta decay experiments (GENIUS, Majorana, 
MOON, EXO) pi] . We finally make a remark that allowing a non- vanishing 5\i will not qualitatively change 
the pattern of our fit. 
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SNO-s-Chooz+SiiperK 




FIG. 3. We present the available region in the 8b = \Sb\ versus 5 M = | (5 M | plane for a combined fit using Chooz, 
SuperK and the SNO results for the MSW-LMA solar oscillation solution, with cos (3 = 1. The circles are solutions for 
6\ ^ 0, the crosses are for 8\ = 0. This figure should be compared with Fig. 5 of ref. [17]. 
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